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Length Scales of Acceleration for Locally Isotropic Turbulence 
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Lengtli scales are determined that govern the behavior at small separations of the correlations of 
fluid-particle acceleration, viscous force, and pressure gradient. The length scales and an associated 
universal constant are quantifled on the basis of published data. The length scale governing pressure 
spectra at high wave numbers is discussed. Fluid-particle acceleration correlation is governed by 
two length scales; one arises from the pressure gradient, the other from the viscous force. 
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Accelerations in turbulent flow are violent ^ and 
are important to many types of studies. Experiments 
in which particles are tracked in turbulence have found 
accelerations as large as 1500 times that of gravity and 
have quantified the probability density of accelerations 
to show that extreme events are likely because the ac- 
celerations are highly intermittent Q. Understand- 
ing turbulent accelerations is essential to many studies, 
includingdispersion and settling of particles in turbu- 
lence other interactions of bubbles and particles 
in turbulence effects on flying insects pollutant 
transport B, and to the applicability of Taylor's hy- 
pothesis fS]^. Turbulent accelerations are the object 
of several studies aimed at solving the mystery of rain 
onset from liquid-water clouds; the accelerations affect 
collision-coalescence, preferential concentration, and lo- 
cal supersaturation ||7||. To advance understanding of 
turbulent accelerations, the length scales of acceleration 
correlations are quantified here, and those scales are dis- 
cussed in connection with the scaling of the pressure spec- 
trum by Gotoh and Fukayama Q and the inertial-range 
scale found for fourth-order velocity structure functions 
by Kerr et al. §. 

Universal scaling of the small-scales of turbulence 
and its analogue in other areas of physics are discussed 
in Ref. The K41 scaling parameters are 

mean energy dissipation rate e and kinematic viscos- 
ity V. Traditional scaling uses K41 scaling of small- 
scale statistics, and seeks the resultant R\ dependence 
of the scaled statistic that Rx dependence is 

akin to symmetry breaking. Here, R\ = (u^^'^'^ \t / v \s 
the Reynolds number, where Taylor's (2) length scale is 
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{u\) / Udujdx^y 
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and (u\^ is the variance 

of one component of velocity. Since (itf) is a prop- 
erty of the large-scale structure, so is i?A- A radical de- 
parture from tradition is to deduce scaling relationships 
from relationships amongst statistics as determined from 
the equations of motion. That method leads, for exam- 
ple, to the assertion that mean-squared pressure-gradient 
scales with an integral of the fourth-order velocity struc- 
ture function for such deduced scaling, there is no 
residual dependence on large-scale parameters Q . The 
deductive method is extended herein to the length scales 



of acceleration. To quantify the deduced length scales 
it is necessary here as in Ref. 



13 



to use existing tur- 
bulence phenomenology; that causes both e and R\ to 
appear herein. However, the deductive theory does not 
explicitly contain K41 scales or Rx. 

The fluid-particle acceleration is related to the 
acceleration caused by the viscous force, vV'^Ui, and that 
caused by the pressure gradient, —dxiP, by the Navier 
Stokes equation: = Dui/Dt — —dx^p + vdx^dx^Ui, 
where Ui is the velocity, p is pressure divided by fluid 
density; density is constant. Here, = dx^dx„ is 
the Laplacian operator, d denotes differentiation with re- 
spect to the subscript variable, summation is implied by 
repeated Roman indices. Consider two points of mea- 
surements, X and x' = X + r; let r = |r|, and let prime 
denote evaluation at x'; e.g., p' = p (x', i), = Ui (x', i), 
etc. Two-point differences are denoted by Ap = p ~ p' , 



Aui 



etc.; angle brackets, i.e. 



denote an 



average. Let the 1-axis be in the direction of r such that 
the 2- and 3-axes are perpendicular to r. Taylor's ||] 
length scale At is the length scale of the parabola that 
osculates the velocity correlation at r = 0. Speciflcally, 
forr = (ri,0,0), (uiO / (u?) = 1 - (ri/2AT)' + • • •. We 
consider the spatial correlation tensors of fluid-particle 
acceleration (0^0^), pressure gradient (^dxiPdx'.p'^ , and 
viscous force per unit mass (^vV'^Uii'V'^iu'j'^ , as well as 
the pressure structure function Dp (r) = (Ap^) and spa- 
tial spectra. We determine the length scales of those 
statistics. The length scales considered in this paper are, 
like At, the length scale of the parabolas that osculate 
those spatial correlations at r = 0. Unlike At, all the 
length scales defined here do not depend on a large-scale 
parameter like (wi). 

The considered statistics, (dxiPdx'.p'^ , Dp, etc., 
are related to the fourth-order velocity structure func- 
tion (i.e., T),jki (r) = (Au^AujAufcAw;)) in Ref. @; 
that theory is based on the Navier Stokes equation, in- 
compressibility, and local isotropy, without further as- 
sumptions. The theory |l^ also relates (i^V^Uii^V^/Wi) 
to the third-order velocity structure function Dj^fe (r) = 
(AuiAujAuk)- The theory has been used to calculate 



dxiPdx'.p') and (^vV'^UivV'^iu'j^ from T>ijki and D 
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respectively, by means of hot-wire anemometry jTsf and 
from direct numerical simulation (DNS) (see Figs. 12, 
13 of Ref. Q), as well as to compare Dp calculated 
from IDijki with Dp calculated from DNS pressure fields 
pTf . The theory also determines the length scales of 
those statistics, which is the topic here. 

Taylor series expansion is used in Ref. to 

formulate the scahng lengths of (dxiPd^'.p'^ ■ For 
brevity, denote the mean-squared pressure gradient by 
X = {dxiPdxiP) (sum on i). The longitudinal component 
of the pressure-gradient correlation is (dxipd^'^p') (the 
1-axis is in the direction of r). Power series expansion 
of the theory's |Q relationship between {dx^pdx'^p') and 
components of D^t;; gives 



{dx^pdx'.p') = (x/3) (1 - rV2A? + ■■■)■ 

nl/2 



x/iQUdx^ui)- 



'[{{dxtU2)} {{dx^Uif {dx^U2) 

1 + x4 r - -, ^ — 
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(cJkUJkSijSij) 
(SijSij) 



{SijSij ) 



(1) 

(2) 

(3) 
(4) 



Thus, Ai is the scaling length for the longitudinal com- 
ponent {dxiPdx^p') , analogous with Taylor's scale. In 
Eq. (|^), hg is given in terms of 3 of the 4 fourth-order 
invariants determined by Siggia |l^ , where Sij is the rate 
of strain tensor and Uk is the vorticity. The invariant 
missing from Eq. ^ is the one that Siggia ||l^ expected 
to vary with Rx differently than the other 3, so there is no 
conflict with the expectation expressed in Ref. that 
Hq becomes a constant at large Reynolds number. Also, 

note that hg = {^(V^p)^^ /1575 (^{stySijfy The cor- 
responding length scale of the transverse component of 
the pressure gradient correlation is 3"'^/^Ai (Ref. iQ Eq. 

(61)), and that of {dxiPdx'.p') is (9/5)^^^ Ai. The analo- 
gous length scale of Dp (r) is v^Ai (Ref. Q Eq.(38), et 
seqq.). The above shows that to cause pressure-gradient 
correlations to coincide near the origin of a graph, one di- 
vides them by x and uses r/Ai on the abscissa; for Dp (r), 
Dp (r) /x^i should be on the ordinate with r/Ai on the 
abscissa. The pressure spectrum is the sine transform 
of the pressure-gradient correlation [Ref. 10, Eqs. (15), 
(22a)]. Because Ai scales the pressure-gradient correla- 
tion at small r, it follows that to cause pressure spectra 
to coincide at viscous-range wave numbers, one should 
divide the spectra by x^i use kXi on the abscissa (fc 
is wave number). A Reynolds number limitation for this 
to be so is discussed near the end. 



To calculate Ai using Eq 
using Eqs. (^ or (|4|). Siggia's [|l9 




we must evaluate hq 
DNS data for Rx - 60 



FIG. 1: The Rx dependence of hq: +, DNS data [gOf; ex- 
perimental data; A, |^|; □, with error bars; solid line, 

iil. 



Kerr pO| and the experimental data of Tsinober et al. 
[21 1, Pearson and Antonia |2^, and Zhou and Antonia 
[23 1, can be used to calculate hq; the result is shown in 
Fig. 1. As suggested by the error bars, the data of Zhou 
and Antonia do not give a reliable value of hq , but those 
data are consistent, within experimental error, with the 
other data. The data of Pearson and Antonia [^ (see 
their Fig. 7), for which error bars are not available, have 
scatter similar to those of Ref. [^. However, by fitting 
their data to straight lines over a broad range of Rx: 
Pearson and Antonia [0 mitigate random error to find 

: 5.52 



{dx^U2f^ / {{dx.Ui)" 



that for 37 < i?A < 10^ 

and Udx.uif {dx,U2f) / (^{dx.uit) = 0.84. Use of 

Eq. (y) then gives hq — 0.32; that value is shown as the 
solid line extending over the range 37 < Rx < lO'^ in Fig. 
1. Note the slight increase of /in with decreasing Rx at 
Rx < 20. The data of Refs. are in such good 

agreement that hq is probably independent of Rx for the 
range 20 < Rx < 10'^, wherein hq ~ 0.3, and this value 
might persist for Rx oo. 

Because Ai is essential for scaling, we next quan- 
tify it by means of existing phenomenology. Use of ex- 
isting phenomenology requires use of the K41 accelera- 
tion scale, e^^'^v'^^^^ , and introduction of r/ = (t/'^/e)^^^, 
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s = Ibv (^{dx^uify and F = ({dx^ui)"^) / ({dx^uif 



into Eq. (g) 
imation, as 



Doing so reexpresses (E^, without approx- 



Ai = (5/2)7? (x/£^/2i.-i/2)'^'(F/jQ) 



-1/2 



(5) 



to 90 gives a single value hq ~ 0.28. The DNS data of 



We emphasize that the theory [Q does not contain K41 
scaling. Use of Eq. (^) leads to asymptotic expressions 
for Al for high and low Rx as follows. On the basis 
of empirical data for the form of Dim {j"), an approxi- 
mation for X is X 'i.lH^£^/'^v-^''^F° '^'^ for Rx > 400 
. This formula produced quantitative agreement [T^ 
with the DNS data of Refs. ||l6| and js) when combined 
with the data for F at high Reynolds numbers from Ref. 
||l|, namely F ~ IMR^^^. Recently, the same DNS 
data as that of Ref. ||] have been used to give val- 
ues of F in Ref. |2|. Those values at Rx > 400 are 
consistent with F ~ l.lSRl^^. To use F ~ lASR^^^ 
instead of F ~ 1.36i?°'^^ and to maintain quantitative 
agreement of x with the DNS data of Refs. [|l6j-[^, it 
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FIG. 2: The dependence of 
Eq. (§); *, DNS data |||| in Eq 
Eq. ®; , Eq. (|). 



rj: +, DNS data 
A, DNS data 



is necessary to use the coefficient 3.5 instead of 3.1 in 
X ^ S.lH^e^^'^iy-'^/'^F" '^^. That is, to be consistent with 
the data of both Refs. [|1 and S, we use, for Rx > 400, 



3/2„-l/2p0.79 



3/2„-l/2^0.25 



(6) 



Evidence that is a constant for Rx > 80 and that its 
value is about 0.65 is given in Ref. [|l3j. For — 0.65, 
Eq. agrees quantitatively with the DNS data in Table 
1 of Ref. § for Rx > 387; also, Eq. (|) is a good 
approximation of the DNS data for Rx as small as 200. 
Substitute Eq. (|) in Eq. (|) to obtain, for Rx > 400, 



Ai ~ 4.777 {HJF' 



0.21; 



.1/2 



6.877i? 



-0.033 



(7) 



= 0.65, 



For the right-most expression in Eq. 
hq — 0.3, and F ~ 1.18i?° '^^ were used. For the case of 
low Reynolds numbers, Ref. |Q obtains 

X O.lle^^^ ly-^/^Rx for Rx < 20, (8) 

which agrees with the DNS data of Ref. For low 

Reynolds numbers, use Eq. (S) in Eq. (0) to obtain 



Ai ~ Qmri[Rx/FhQf''^ for Rx < 20 



(9) 



None of Eqs. to is K41 scaling because of their 
Rx dependence. Figure 2 shows Ai increasing relative 
to rj as Rx increases from 9 to attain a maximum near 
Rx ~ 150, beyond which Ai gradually decreases relative 
to rj in agreement with Eq. (|7|). 

We now consider the spatial correlation of the vis- 
cous force because it it is part of the fluid-particle cor- 
relation as well as of intrinsic interest jlj]. The spatial 
correlation of z^Vju is l^alial: 



(r) ^ (i.2v2«,v2,«0 = _^V?V?A, (r) (10) 



W4)V2a,.,A,/c (r). 



(11) 



where is the Laplacian operator in r-space. Then, 
for the longitudinal components of Vij (r) , Eq. ( pd] ) yields 

Viiir) = (z//2r3) [i?ni (r) + 2i?i22 (r) 

-5rdrDi22 (r) ~ r^dlD^22 (r) ] . (12a) 



A similar formula applies to the transverse component 
^22 {r) jl^. Power series expansion gives Dm (r) = 

{d^,u,f) [1 - (rV (2A?,^^ J) -f • • •] , where 



-4((9,,ui)' 



Then, Eq. ( tL2a| ) yields 
35 
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(13) 



^11 {r) = -y^ {{d.^u.f) [1 - (rV2A^,J -f 



where Ay^^ = A£)^^^/\/18/5 is the length scale of the os- 
culating parabola of Vn (r). The corresponding expan- 
sions of the transverse component of Vij (r) and Vu (r) 
give length scales Xv^^/V^ and Xvn/ \/^Ji, respectively. 
Kolmogorov's equation can relate the velocity-derivative 
moments in Eq. ( [T^ ) to other high-order derivative mo- 
ments; this gives the same result as would beginning with 



the relationship of Vij (r) to Dij (r) in Eq. (10). The 
difhculty of observing the initial fall-off of Vu (r), and 
hence its length scale Av\i is clearly illustrated in Fig. 3 
of Ref. ^ and Fig. 13 of Ref. ||l|| and their discus- 
sions; one requires a spatial resolution substantially finer 
than rj. Such fine resolution data is not yet available. 
The high-order derivative moments in the denominator 
of Eq. (^3|) likewise require fine spatial resolution; their 
Reynolds-number dependence has not been investigated, 
but it would be surprising if they obey K41 scaling. 

Despite these difficulties, we can obtain the r- 
values at which Vu (r) /Vu (0) = 0.5, and similarly for 
the transverse component, from the graphs of those cor- 
relation coefficients in Fig. 13 of Ref. [1^; the r- values 
are 3.5?? and 2.577 for longitudinal and transverse compo- 
nents, respectively. Limited to Rx = 230 of the DNS 
run in Fig. 13 of Ref. jl^, 3.577 and 2.577 are the first 
estimates of these length scales. It is not suggested that 
these length scales maintain fixed ratio to rj as Rx varies. 

The fluid-particle acceleration correlation is 



A, 



(r) EE (a,a;.) = id^^pd,,p') + V,, (r) 



(14) 



As a result of the differing length scales of (^dxiPd^'P 

and Vij (r) , the theory requires that the initial fall-off of 
the longitudinal component An (r) be described by two 
length scales; the same is true for the transverse com- 
ponent and for An (r). For example, consider the trace 
An (r) and Rx — 230. From the above discussion we 
have Xvii — 3.577 such that Xy^ = Xvn/ ^/5/3 = 2.777. 
From Fig. 2 at Rx = 230 we have that the length scale of 
{d,,pd,>p') is v^Ai = 7.677. Also, Vu (0) /x = 0.015 
at Rx = 230 (TsUll. Thus, A^i{r) has an initial 
rapid but small-amplitude decay (1.5% of the total) from 
Vii (r), having scale 2.7?7, followed by the larger ampli- 
tude and more gradual decay of {dx^pdx'p') with scale 
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7.6r]. It is not implied that the two length scales main- 
tain fixed ratio as R\ varies; the example above is for 
Rx = 230. 

We noted above that to cause pressure spectra to 
coincide at viscous-range wave numbers, one should di- 
vide the spectra by x^i ^^^id use fcAi on the abscissa. 
However, the Reynolds number must be large enough 
that, at high wave numbers, the pressure spectrum has 
negligible contributions from the sine transform of the 
pressure-gradient correlation at r on the order of the in- 
tegral scale. How large is large enough? From the 
pressure spectra from isotropic DNS shown in Fig. 5 of 
Gotoh and Fukayama , where krj is the abscissa, com- 
bined with the weak variation of Xi/rj shown in Fig. 2 
for 100 < i?A < 400, it seems that Rx > 200 is large 
enough. Of course, no such limitation applies to scaling 

of (dxiPdx'p'^ and Dp (r) with the parameters x ^-nd Ai. 

Kerr et al. Q identify, by empirical means, a 
length scale within the inertial range of fourth-order 
velocity structure functions {T>ijki (r)) such that scal- 
ing exponents should be determined only from r greater 
than that scale. In the inertial range, the divergence 
of Dyfci (r) equals the pressure-gradient velocity-velocity 



structure function which is, in turn, related to Aij (r) 
(Eqs. 9, 10, A4, A5 of Ref. ||5|). Taylor series expan- 
sion of those relationships shows that Ai is a scale of, at 
least, linear combinations of the nonzero components of 
Dyfei (r). From Fig. 2 of Ref. their empirical scale 
is about 5 times Ai; on the other hand, Ai is only the 

initial roll-off of (dxiPdx'.p'^ ■ Further work is needed to 
establish a causal relationship between the scales, if one 
exists. 

The theory jlj] predicts the scaling with x that 
was found empirically in Refs. ||l^ and Even for 

large Reynolds numbers, dividingl^dxipdx'^p') by x ^nd 
Dp (r) by x-^i and similarly for pressure spectra is not 
K41 scaling. The approximation Eq. contains, in 
addition to K41 scaling parameters, the factor F and 
thereby dependence on Rx. The deviation of both x 
and F from K41 scaling is large for large variations of 
Reynolds number. Whether or not Xi/rj continues the 
weak downward trend shown in Fig. 2 or becomes con- 
stant as Rx ^ oo is also a topic for further investigation, 
as is the relationship of r] to the viscous-force scale Xvn ■ 
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